Quick note

Have transfered Grodent's atmospheres, total H2 ionisation rates and subsequent H3+ densities into a spreadsheet. The first bit needs conversion so the five species can be used by RIDE properly. The second bit is the direct comparison between my ionisation rates and his. The final bit compares our chemistry models. Might be interesting...

I was a little worried that leakage of the 20keV maxwellian beyond RIDE's upper limit would cause problems with electrons not creating enough ionisation at the bottom of the atmosphere. However, the atmosphere cuts off before that point anyway...

More and more and more...

Another block of maths demonstrating out of the way.

One run has finished. Saturn1000 slot 100keV GS1m run. Analysis was as expected. Recoded as GS10 version of same energy and set off again. Have decided to kill off all the fast runs asap rather than do them concurrently with the slow runs, then do all the slow runs in one big block - asides from those currently being done and their heirs...

Not that any 'fast' runs are left asides fro the 100keV ones. All medium to slow.

Also had a quick look at Makenzie's newly produced profiles. She hasn't produced the further out ones, which are the important ones yet, so will wait until tomorrow for discussions. Or later. Who knows...

...one last thing. References for Martian cross-sections kept from Europlanet, this list I will be going through for obvious reasons:
Straub et al 1996 - already got
Lindsay et al 2000 - now got
Manyan et al 2000 - ah, ManGan et al... ahem, now got...
Deutsch et al 2002 - now got
Haider et al 1992 (modelling) - JGR, dammit!
Erdman and Zopt 1983 (cameron band?) - ZIPF, what sort of a name is that? Far prefer Zopt... Anyway, not included in UCL's subscription to Planetary and space science, so no-go despite importance...
Fox and Dalgarno 1979 - JGR - though did get hold of three 1979 papers by Dalgarno relating to proton processes I was in dire need of for DICE...
Lummerzhien and Lilensten 1994 - Lummerzheim, given I own half his back catalogue, it sometimes seems, I should know how to spell his name... plus annales has withdrawn the PDF! Did they know I wanted it?
Tran and Vidal 1998 - TIan... now got, oh!, there's another, got it too. Waaaaait. A third useful Tian and Vidal 1998 has appeared with cross-sections, my they had a busy year - nope, here's number four, maybe Tian and Vidal 1998 is a generic reference to their entire career...
Torr and Torr 1985 (Photon)
Avalgan et al 1998 (Photon)
Hitchcock et al 1980 (Photon)
Masuoka 1994 (Photon)

Will deal with the photon stuff another time. Mostly have them, but no real prospect of a photon thing to deal with yet as Tommi is otherwise engaged.

I make no garrentees that any of those names are spelt correctly, or in any way resemble the correct spellings.

All qualified stuff today

ie, looking at programs, program cycles etc and deciding what to do.

RIDE produces heating rates and ion production rates. The effect on this on the atmosphere needs to be quantified. The initial problem is the H3+ column densities supported by auroral events. This means taking the production rates and sticking them into chemical rate programs before iterating through to produce H3+ densities, then integrating along the line of sight. RIDE produces production rates, chemical reaction rates module has come from KIM, and can be used with improvements from JIM and other sources. But this is just stage one...

The heating through chemistry and auroral events needs to be transported, as does the new set of chemicals. The heating also needs to iterate back into the rates. One part of the greater production rate at the peak of the ion production profile is that there is also a close by peak of the heating rates, allowing a greater number of endothermic reactions to take place.

When did I get into atmospheric chemistry? Or physics for that matter...?

To do this, need to understand the chemical reaction rates routine... make it my own... hehe

Inputs:

dtimeplas,iwrite,in,is, - related to position in flux tube, need to be rewritten for altitudes
He,H,H2,CH4, - neutral densities from RIDE, need to be expanded, eventually
Hp,H3p,Hep,elec, - ionic densities from RIDE - as above
prodhp,prodh3p,prodhep, - production rates of certain things from ride
betaChp,betaCh3p,betaChep,betaRhp,betaRh3p,betaRhep, - reaction rates, initial
CRC,vib_h2_sw,vib_h2_param, - related parameters
tarr,rch3p, - recombination rates
temp,recomb,chex - temps, plus charge exchange rate?

Usual declarations and inclusions. Most the things to alter from flux tubes to altitudes here.

Check for presence of specific things

All rates are zeroed to begin with

Rates are then calculated as a power law. Not all rates included, but pointers to remaining JIM rates are given where appropriate, and obviously they set up a pattern that can then be used for other rates beyond the original scope of this thing.

Rates then stuck into tubes and corrections for underflows etc made.

Have also been digging around for other things such as photochemistry et al, which provide either useful tools or contain bits of code that make the chemistry more complete - eg changing production and loss rates into actual chemical changes.

Opened thesis, didn't edit it though...